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Abstract 

In this short review we present our recent results concerning the rotation of 
atomic Bose-Einstein condensates confined in quadratic or quartic potentials, and 
give an overview of the field. We first describe the procedure used to set an atomic 
gas in rotation and briefly discuss the physics of condensates containing a single 
vortex line. We then address the regime of fast rotation in harmonic traps, where the 
rotation frequency is close to the trapping frequency. In this limit the Landau Level 
formalism is well suited to describe the system. The problem of the condensation 
temperature of a fast rotating gas is discussed, as well as the equilibrium shape of 
the cloud and the structure of the vortex lattice. Finally we review results obtained 
with a quadratic + quartic potential, which allows to study a regime where the 
rotation frequency is equal to or larger than the harmonic trapping frequency. 

The possibility to obtain quantum degenerate gases by a combination of laser and 
evaporative cooling has opened several new lines of research, at the border of atomic, 
statistical and condensed matter physics (for a review, see e.g. [H El El Hj ) • Among them, 
the rotation of a Bose-Einstein condensate raises many interesting problems with respect 
to the case of a classical fluid. Since the condensate is described by a macroscopic wave 
function tp(r) = p{r)e l ^ r \ where p and are the spatial density and phase of the fluid, 
there exist strong constraints on the velocity field of the rotating gas. In a place where the 
spatial density is not zero, this velocity field is given by v = h'V^/M (M is the particle 
mass), hence V x v = 0. The circulation of the velocity field is quantized along any close 
contour on which p ^ 0, and it is a multiple of h/M. The rotation of the fluid is thus 
only possible through the nucleation of quantized vortices [2J E] , which are singular points 
(in 2 dimensions) or lines (in 3 dimensions) of vanishing density, and around which the 
phase varies by multiples of 2tt. Vortices are universal objects which appear in many 
macroscopic quantum systems, such as superconductors and superfluid liquid helium. 

"Unite de Recherche de l'Ecole normale superieure et de l'Universite Pierre et Marie Curie, associee 
au CNRS. 
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Since the achievement of Bose-Einstein condensation in atomic gases, many experi- 
mental and theoretical studies have been devoted to vortices in these systems ([7] and refs. 
therein). A typical experiment is the following: One starts with a condensate initially 
at rest, confined in an axisymmetric trap (symmetry axis z), and stirs it by applying an 
elliptic potential rotating at frequency Q around z. For very small values of Q no angu- 
lar momentum is transferred to the condensate. Just above a critical value Q c , a single 
vortex is nucleated jSj. For a system in thermal equilibrium, the existence of this critical 
frequency can be viewed as a manifestation of superfluidity: for a slow enough rotation 
frequency, the stirrer cannot drag the condensate and set it in motion. For stirring fre- 
quencies notably larger than Q c , the number of vortices in the condensate N v increases, 
and values up to N v = 200 have been obtained experimentally P ITUj. 

In the following we will be mostly interested in the large vortex number case, which 
is achieved by choosing Q close to the trapping frequency u± in the xy plane. Note that 
in a purely harmonic trap, one must keep Q below uj±] the centrifugal force otherwise 
exceeds the trapping force and the gas is destabilized [TT]. The main features of the 
vortex assembly in this large N v regime are well known. The vortices form a triangular 
Abrikosov lattice with a surface density n v = MQ/(irh) [T2]. When Q — > u±, the radius of 
the gas tends to infinity since the confinement by the trapping potential is nearly balanced 
by the centrifugal force. Since the surface density of vortices is constant (~ Muj±/(nh)), 
the number of vortices N v also increases to arbitrarily large values. 

In principle the number of vortices can reach and even go beyond the atom number 
N. For such a fast rotation, the description of the system by a single macroscopic wave 
function is expected to fail, and the ground state of the system should be strongly corre- 
lated. Up to now this regime has been investigated only theoretically [T3J EH EES UH UZ] 
and will not be addressed here. 

The outline of this short review is as follows. In section^ we present the mechanism 
used to set a condensate in rotation. Then in section|2]we briefly review some experiments 
performed with a single vortex condensate. In section |SJ we focus on the fast rotation 
regime and discuss some important features of this system, such as its condensation 
temperature and its equilibrium shape. We use the Landau Level approach, which makes 
a nice connection between this physical problem and that of charged particles in a uniform 
magnetic field. In section 0] we turn to a configuration that we recently investigated in our 
laboratory, which consists in superimposing a trapping quartic potential onto the usual 
quadratic one. This allows us to explore the rotation regime Q > u>j_, and we review some 
results obtained for the vortex patterns in this regime. We conclude in section El by giving 
some perspectives of this rapidly evolving field of research. 

Note that due to the lack of space this paper does not attempt to be a full review of 
the work that has been performed on rotating quantum gases but will be subject to the 
following restrictions: First we will only discuss studies on single component condensates 
and refer the reader interested in rotating spinor condensates to the recent experiments of 
the Boulder group ( |18j and references therein). Second we will illustrate our discussion 
using mainly experimental results from our laboratory and only give references to the 
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achievements of other groups. Finally our reference list focuses on papers published after 
2001. For a discussion of earlier articles see the detailed review paper by Fetter and 
Svidzinski [7j. 



1 Setting a BEC in rotation 



In order to nucleate vortices in a condensate, two classes of methods have been used. The 
first one consists in imprinting on the condensate the phase pattern e l6 of the desired 
wave function. It was successfully implemented experimentally by the Boulder group for 
a two-component condensate jTHj. A related scheme, based on the adiabatic inversion of 
the magnetic field at the center of the magnetic trap, has been used at MIT The 
second method, which is used in our group, consists in using a mechanical stirring of the 
condensate. In order to do this, one can use the potential created by a laser |HJ [§] or by a 
magnetic field |2H I22|. While the phase-imprinting method is well suited for nucleating 
a single vortex, the stirring approach seems to be more flexible and allows nucleation of 
a large number of vortices. 

Once the gas is rotating, a third method can be used to increase the angular momentum 
per particle. It consists in eliminating atoms with an angular momentum smaller than the 
average, so that the remaining particles rotate at a larger angular speed. This "evaporative 
spinup" method has been implemented in Boulder |2HJ. 

We now present the system that we have been using for rotation experiments in our 
group. We use rubidium ( 87 Rb) Bose-Einstein condensates produced in a cylindrically 
symmetric Ioffe-Pritchard trap, with a frequency u± in the xy plane and u z along the z 
axis. The magnetic trapping potential thus reads: 



Typically oo± ~ 10 u z in our experiments so that the equilibrium shape of the condensate 
is an elongated cigar. For uj z /2tt ~ 10 Hz and iV~3x 10 5 rubidium atoms in the trap, 
the length of the cigar is 100 /xm and its diameter is 10 fim. 

We stir the condensate with a laser beam propagating along the z axis. The beam has 
an anisotropic cross section and its eigenaxes rotate at a frequency Q. The time-dependent 
potential created by the laser beam can be written as 



The parameter e is a dimensionless measure of the relative strength of the stirring and 
the magnetic potentials. In practice we choose e ~ 2-10 %. 

We apply the stirring potential onto the condensate for a fraction of a second, in 
order to transfer angular momentum to the gas. The condensate then equilibrates in the 
cylindrically symmetric potential (0) for ~ 1 second. The trapping magnetic field is then 





V stiT (t) = e -Mu\ [(. 



x 2 - y 2 ) cos(2fit) + 2xysm(2ttt)] 
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Figure 1: Quantized vortices. Absorption images of a stirred Rb Bose-Einstein condensate. 
The rotation frequency is increasing from left to right (for details see jH]). 

switched off and the gas undergoes a ballistic expansion for a period of ~ 20 ms. Finally, 
we perform absorption imaging along the rotation axis z. The vortices which have been 
nucleated in this process appear as density dips in the images, as seen in FigJTJ 

When studying theoretically the problem of a rotating gas, one usually assumes that 
an arbitrarily small anisotropic potential, rotating at frequency Q, is added to the main 
isotropic trapping potential. In presence of this stirring potential, the frame rotating at Q 
is the only one in which the state of the system is stationary. The hamiltonian H in this 
rotating frame is deduced from the hamiltonian in the lab frame by H = — QL Z . 
Experimentally, as we just described, the rotating potential is often switched off for some 
period before the measurement. Q then plays the role of the Lagrange multiplier associated 
with the deposited angular momentum L z , which is a constant of motion when the system 
evolves in the axi-symmetric potential. 

In the following we will mainly focus on the equilibrium properties of the rotating 
system. We refer the reader interested in the dynamics of vortex nucleation and decay to 

I2H ESI [2E1 EH [2H1 [2H E01 EU E21 E31 EH ESI EEl E3 EH 
2 Single vortex physics 

For a proper choice of Q and the equilibration time after stirring, it is possible to nucleate 
in a reliable way a single vortex in the center of the condensate. Several experimental 
studies have been performed on such a system. First, the average angular momentum 
per particle L z has been measured and found to be of the order of Ti [39J. This measure- 
ment was performed using the relation between L z and the frequencies of the transverse 
quadrupole modes of the condensate jlU] • Complementary information has been obtained 
using atom interferometry to measure the phase pattern of the wave function [4*T| I42j . 

Concerning the vortex line itself, its equilibrium shape has been determined: the line 
is often curved at the two ends of the cigar [33]. This bending is a symmetry breaking 
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Figure 2: Kelvin mode of a vortex line. Transverse images of a Bose-Einstein conden- 
sate with a single, positively charged vortex, obtained after excitation of the transverse 
quadrupole mode m = —2 (a) and m = +2 (b). Fig. |2K gives an evidence for the Kelvin 
mode of the vortex line. This mode has an angular momentum m — — 1, and it is thus ex- 
cited by the decay of the transverse quadrupole mode m = —2 into two "kelvons" (quanta 
of the Kelvin mode). The two kelvons have the same energy and propagate in opposite 
directions, which ensures the conservation of linear momentum. By contrast, Fig. EJd 
shows no oscillation of the vortex line, as expected since the decay of the quadrupole 
mode m = +2 into kelvons (m = — 1) is forbidden by angular momentum conservation. 
Figs. Etc,d) are the corresponding density profiles (for details see |Hlj). 

effect, and it can be understood by noticing that a radially centered vortex is favored 
at the axial center of the cigar, where the density is large, whereas it costs less energy 
to radially off-center the vortex line at the ends of the cigar, where the density is low 
|H| 0H1 113 HH EH] • Some specific normal modes of the vortex line have also been 
observed, such as the precession of a single vortex when it is not aligned with the trap 
axes EHJEHI521I53I, and the Kelvin mode of the vortex line (see Fig. |2D [Ml ESI EH1 E3 EE| • 

3 Fast rotation in a harmonic potential 

We now address the case of a fast rotating gas where the average angular momentum per 
particle is large compared to H, i.e. in which many vortices have been nucleated. We 
consider that the gas is confined in a purely harmonic potential We first present the 
Landau level approach to this problem, which is directly connected to the description of 
the motion of a charged particle in a uniform magnetic field [53 EH1 HH] • We then address 
the determination of the critical temperature of the rotating gas, and we compare the 
result of the semi-classical approach given in [62] with the treatment using the Landau 
level basis. We then turn to the discussion of the equilibrium shape of the rotating 
condensate and the structure of the vortex pattern. 
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Figure 3: Landau level structure. Single particle spectrum for the transverse motion (in 
the xy plane) for Q = 0.9. The index k labels the Landau levels. The energy is expressed 
in units of hu±. 

3.1 The Landau level approach 



In the frame rotating at frequency Q, the non-interacting, single particle hamiltonian is 
H = H ± + H z with H z = P 2 /(2M) + Mu 2 z 2 /2 and 



p2 + p2 1 

-^^ + ^Mu 2 ± (x 2 + y 2 )-QL z 



2M 



(3) 
(4) 



where the vector potential is A = Mfl x r. Eq. (jlj) is formally identical to the hamiltonian 
of a particle of charge 1 placed in a uniform magnetic field 2MQz, and confined in a 
harmonic potential of frequency \Joj\ — Vt 2 . The Coriolis force, which has the same 
mathematical structure as the Lorentz force, originates from the vector potential A, 
whereas the term — MQ 2 (x 2 + y 2 )/2 corresponds to the centrifugal potential. 



Common eigenstates of H and L z have single particle energies 



Ej,k,n/h = oj± + — + j(u_L — Q) + k(u± + fl) + nu z 



(5) 



and angular momentum h(j — k), where j,k,n are non-negative integers. For Q close 
to uj±, the contribution of the transverse motion to these energy levels (terms in j and 
k) groups in series of states with a given k, corresponding to the well known Landau 
levels (Fig. EJ). The lowest energy states of two adjacent Landau levels are separated by 
h(u±_ + Q), whereas the distance between two adjacent states in a given Landau level is 
h(ujj_ — Q). When Q = u±, all states in a given Landau level are degenerate. Physically, 
this corresponds to the case where the centrifugal potential exactly balances the trapping 
force in the xy plane, and only the Coriolis force remains. The system is thus invariant 
under translation, hence the macroscopic degeneracy. 

When interactions between particles are taken into account, the Landau levels are no 
longer eigenstates of the N— body hamiltonian. However they are still relevant in the 
regime of fast rotation. Indeed as Q — > ujj_, the restoring force in the xy plane becomes 
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very small and the density of the gas drops. The interaction energy per particle is then 
small compared to the distance 2hu± between two Landau levels, and the states of interest 
are essentially those associated with k — 0, i.e. the lowest Landau level (LLL) [qTJIIqTH IHT] . 
Any function ip(x, y) of the LLL can be cast in the form: 

i>{x, y) = e -(- 2 +y 2 )/^l p( x + i y ) , (6) 

where a\ = h/(muj±) and P(u) is a polynomial (or other analytic function) of u. When 
P(u) is a polynomial of degree n, it has n complex zeroes. Each zero is the position of 
a singly-charged, positive vortex, since the phase of if)(r) changes by 2ir along a closed 
contour encircling the zero. 



3.2 Critical temperature for a rotating gas 

The critical temperature for an ideal rotating gas in a harmonic potential has been derived 
by S. Stringari using a semi-classical approach jS2]. Let us briefly outline the reasoning. 
One starts from the semi-classical relation between the atom number N and the temper- 
ature T at the BEC transition: 

N = J_ f £ r d s \ , (7) 

h A J exp(H (r,p)/k B T) - 1 

where we assume that the minimum of the trapping potential is at zero energy, so that 
we set the chemical potential also equal to zero at the transition point. Using the form 
(JU of H± and making the change of variables p' x = p x + MQy, p' y = p y — MQx, we obtain 
a new integral. This integral is identical to the one giving the condensation criterion for 
a gas at rest, confined in a cylindrically symmetric, harmonic potential. The transverse 
and longitudinal frequencies for this model system are \Ju\ — fi 2 and io z , respectively, 
and we thus get: 

kBT\ 3 -.3/2 



N = c{3) \^) ' ^ = K-n>-, (8) 

where ((x) = Yl n n ~ x ^ s ^ e Riemann zeta function (C(3) — 1.202). This entails that, 
within the semi-classical approximation, the Coriolis force associated with the vector 
potential A(r) has no effect on the critical temperature. This is formally identical to the 
Bohr - van Leeuven theorem, stating that there is no magnetism at thermal equilibrium 
in a system of charges described by classical mechanics. The only effect of rotation in (jH| 
is the change of the transverse frequency due to the centrifugal force. 

The result (JH| can be recovered from the exact one-body spectrum in terms of Landau 
levels given in ( J^j l. Indeed this spectrum is the same as that of a 3D harmonic oscillator 
with frequencies u± — Q, u± + Q, u z , whose geometrical mean is the frequency u, hence 
the result ©. 

We now discuss briefly the validity condition of (jH|. From the reasoning based on the 
Landau level structure, we see that k B T must be large compared to each of the three 
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energies h(u± ± Q) and tku z . In practice, we have uj±_ ^> c<j 2 for a cigar-shape trap, so that 
when Q ^ the most stringent validity condition for the use of (JBJ) is ksT ^> 2fru;_i_. 
It is not easy to recover this validity condition directly from the use of (jHJ). Indeed one 
could have naively expected that k B T ^> h^Juj 2 ^ — f2 2 , %uj z would be a sufficient condition, 
which is clearly not the case. 

3.3 The equilibrium shape of the rotating condensate 

We now suppose that the gas is at zero temperature and review possible approaches for 
determining its ground state when repulsive interactions are taken into account. As usual 
in the physics of cold gases, these interactions are assumed to be point-like, and they 
are characterized by the s-wave scattering length a s . The ground state 0(r) of the gas is 
obtained in the mean-field approximation by minimizing the energy per particle [3J H] 

E[<f>} = J (V W] + ^|0| 4 ) d 3 r (9) 

where g = 4irh 2 a s /M and is normalized to unity. 

3.3.1 Rotational hydrodynamics approach 

This approach is based on the approximation of diffused vorticity, where the singularities 
of the velocity field v(r) and of the atom density n(r) at each vortex core are averaged 
out. This method is adequate for describing the system at macroscopic distances, larger 
than the intervortex spacing. From the equations of motion of rotational hydrodynamics 
EH! ( see a ^ so one derives the steady state velocity field v(r) = f2 x r and the 

spatial density n(r) = N\(j)(r)\ 2 : 

n{r) = i max f , [i - V^gO) + ^-{x 2 + y 2 ) \ (10) 

where \x is the chemical potential. This density profile is the usual inverted parabola 
corresponding to the Thomas- Fermi result, for an axisymmetric potential with frequencies 
\Jlo 2 ± — Vt 2 in the xy— plane and uo z along the z axis. 

The result (HOj) is valid only if fi ^> Two z . In the opposite regime \x <C hu: z , the z motion 
is "frozen" to its ground state (a gaussian of extension a z = ^JUj (Mu z ) ). The relevant 
wavefunctions can be written as 4>{r) = ip(x,y) e~ z ' 1 l 2a ^ and one has to minimize 

EM = J U* [H ± 1>] + ^l^l 4 ) d 2 r (11) 

where G = g/(V2n a z ). The spatial density in the transverse plane n(x,y) = N\vp(x, y)\ 2 
is then: 

n(x,y) = ^max( 0, -\m(u 2 l - tt 2 )(x 2 + y 2 ) J , (12) 



S 



As for the determination of the critical temperature, only the centrifugal potential is 
important in this approximation. The Coriolis force plays no role in the global equilibrium 
shape of the condensate. 



3.3.2 Equilibrium shape in the LLL 

We now suppose that the interaction strength is small enough so that the ground state 
of the system is essentially a LLL wave function, corresponding to the quantum number 
k = in Eq. 03 (/i 2hu±). We also assume that fi <C tku z so that the z motion is 
"frozen" to its ground state (n = in Eq. The use of LLL wave functions allows 

to notably simplify the energy functional in Eq. (jllj) . One can indeed prove after some 
algebra the two equalities: 

(£kin> = (E ho ) = ^ + ^ / r [L z ip] d 2 r (13) 



2 2 

where the kinetic and harmonic oscillator energies are: 

(£m»> = m\ 2 d\ (E ho ) = J r 2 |V>| 2 d 2 r . (14) 

The energy is then given by 



E[ij}]=m + J \Muj ± (uj ± -D,)r 2 \'ip\ 2 + 



NG, 



4 j2 



drr . (15) 



We can express the distances and the energies in units of a± = ^h/Muj^ and huj±, 
respectively. We then find that the minimization depends only on the dimensionless 
parameter jM] 

, T MG uj\ i — , T a s uji 
A = N — = ^— = V8tt N 



fi u;± — Q a z uo^ — Vt 

When A < 1, the interaction term AG|^| 4 plays a negligible role and the minimizing 
function is essentially the ground state of the one-body hamiltonian j = k = n = 0. For 
A 3> 1 the minimum energy state is a linear combination of several states corresponding 
to different quantum numbers j's, and it involves several vortices in the region where the 
atomic density is significant. 



The minimization of (|15p within the LLL has recently been discussed in pB| E3 
Let us briefly sketch the main results. One first defines the coarse-grain average n(x, y) 
of the spatial density n(x,y) = N\ip(x, y)\ 2 , in order to smooth the rapid variations at 
the vortex cores. The energy functional (|T3jl can be written in terms of n instead of 
n, provided the interaction parameter G is renormalized to bG, where b ~ 1.16 is the 
so-called Abrikosov parameter [68J. This parameter arises from the discreteness of the 
vortex distribution: since the wave function ip(x,y) must vanish at the vortex location, 
the average value of |^| 4 over the unit cell, hence the interaction energy, is larger than 
the result obtained if was quasi-uniform over the cell. Once this renormalization of G 
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has been performed, the minimization can be performed by letting n vary over the whole 
space of normalisable functions, the only constraint being that n varies smoothly over a±, 
which is the characteristic length scale for vortex spacing. One finds that the coarse-grain 
average of the spatial distribution is the inverted parabola : 

n(x, y) = — max ( , n - Mu ± (u± - O) {x 2 + y 2 ) ) (16) 

This result is valid when the chemical potential \i is much smaller than the distance 2hu± 
between the LLL and the first excited LL, which amounts to: 



Q-z u>± — Q 

Except for the Abrikosov coefficient b, the two results (|12jl and ([TH]) nearly coincide in the 
fast rotation limit, since io\ — Q 2 ~ 2ujj_(uj]_ — Q) when Q ~ u±. 

The fact that the equilibrium shape of the condensate remains an inverted parabola 
even when the dynamics is restricted to the LLL has been checked experimentally by the 
Boulder group [SH] HO] • 



3.3.3 Structure of the vortex pattern 

In first approximation the surface density n v of vortices in a fast rotating condensate is 
uniform. One can show in this case that the coarse-grain average of the velocity field 
is equal to the rigid-body rotation result v(r) = Qz x r ^2], with Q = nhn v /M. For 
rotating BECs, this relation has been checked experimentally at MIT [23]. The vortices 
form a triangular lattice which is known to minimize the interaction energy g J \ip\ A |72| . 

A closer analysis of the vortex distribution shows that the vortex distribution is dis- 
torted on the edges of the condensate [7^1 ESI EZ1 EH] ■ The distortion is particularly clear 
in the LLL, as it can be seen in Fig. 0] obtained by Aftalion et al. j64j . where an example 
of vortex distribution is given for the particular case A = 3000. This distortion of the 
vortex lattice is essential to ensure the proper decay of the atomic density given in ([160. 
Indeed an LLL wave function with a uniform vortex lattice always leads to a Gaussian 
average distribution n(x,y) [61j, instead of the predicted and observed inverted parabola 

(USD. 

Another interesting characteristic of the LLL is that the vortex core is of the same 
size as the distance between adjacent vortices (~ a±). In this respect the entrance in the 
LLL for a rotating condensate in a magnetic trap is the equivalent of the field H C 2 in a 
type II superconductor [71] I65j. Finally we note that the dynamics of the vortex lattice 
itself raises many interesting problems. In particular the so-called Tkachenko modes of 
the lattice have been analyzed theoretically [7^] [7_7J [7B1 [7JH EDI EH E2j and observed 
experimentally [7Hj . 
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Figure 4: Vortices in the LLL. Equilibrium vortex pattern obtained by minimizing the 
energy within the LLL, for A = 3000 (figure extracted from jHll)- The distances are 
expressed in units of a±. The LLL trial wave functions have 52 vortices and the circle 
represents the border of the Thomas- Fermi distribution (|16|). 

4 Fast rotation in a quadratic+quartic potential 

In this section we discuss some results obtained with a quartic potential 7r 4 _ added to 
the usual harmonic confinement (we set r\ = x 2 + y 2 ). This quartic confinement allows 
to study the regime of rotation where Q > u±. This regime is unreachable otherwise, 
since the expelling centrifugal potential — MVL^r 2 ^ would exceed the trapping potential. 
The properties of the rotating gas in a quartic potential have recently attracted a lot of 
theoretical attention |HSl IHU ESI EHl EZl IHH1 IHH1 EDI 

4.1 Implementation of a quartic potential 

Experimentally, we have created a quartic potential using a far detuned laser beam 
(wavelength 532 nm), propagating along the axis of the trap [51] • The waist w of the 
beam is larger that the condensate radius, so that the potential created by the laser 
Uq exp(— 2r 2 L /w 2 ) can be written as: 



The laser frequency is larger than the atom resonance frequency, so Uq > 0. The second 
term in (J 17)) leads to a reduction of the transverse trapping frequency u± and the third 
term provides the desired quartic confinement, with 7 = 2Uq/w 4 . In our experiments 
7 = 6.5 x 10™ 12 Jm" 4 . 



U(r) ~ [/( 




(17) 
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The total trapping potential in the xy plane can be written as: 

V(r ± ) = (\ + 7 , (18) 

where the dimensionless number 7 = 2h , y/(M 2 uj^_) characterizes the relative strength of 
the quartic and the quadratic potentials. For our setup we have u±/2tt = 65 Hz and 
7 ~ 10~ 3 . Hence the quartic term is only a small perturbation of the ground state of 
the one-body hamiltonian. Of course its importance grows when one considers large L z 
states, in which the particle is localized further away than a± from the center of the trap. 



4.2 Critical temperature 

We now determine the critical temperature T c for a gas rotating in a quadratic+quartic 
potential. For simplicity we consider the case Q = u± and we use the semiclassical 
approximation, which is valid if ksT c ^> huj_. Inserting 

Ho(r,p) = ^ + \muj\z 2 + jri (19) 

in (J7J) we obtain 

(20, 

with C(5/2) ~ 1.342. The experimental results shown here were obtained with iV = 3x 10 5 
rubidium atoms. This corresponds to a critical temperature of T c = 60 nK for Q = ou±, 
to be compared with T c = 120 nK for a non-rotating gas. 

Our experiments were performed in presence of radio-frequency evaporation, which 
removes all atoms at a distance r_|_ larger than x ev = 19 /im from the center (this corre- 
sponds to an angular momentum value m = xl v /a 2 ± ~ 200). For Q = uj± the well depth is 
thus Uo = 7xf v ~ 60 nK, similar to ksT c . Since the effective temperature T in evaporative 
cooling is a small fraction of the well depth (typically Uo/ksT ~ 5-10), the rotating gas 
is clearly in the degenerate regime when Q = u±. 



4.3 Observed vortex patterns 

We show in Fig. EJthe images of the rotating gas as the stirring frequency Q is increased. 
For Q < uj±, the vortex lattice is clearly visible. However when Q > u± the visibility of 
the vortices decreases and nearly vanishes for Q = 1.05 uj±_ (= 2ir x 68 Hz). 

The most plausible explanation of this effect is that the vortex lines are still present, 
but strongly bent when Q > u±. This bending may occur because of the finite temperature 
of the gas. A recent theoretical study [88J seems to favor this hypothesis: when looking 
for the ground state of the system using imaginary time evolution of the Gross-Pitaevskii 
equation, it was found that much longer imaginary times were required to reach a well 
ordered vortex lattice for Q > u± than for Q < u±. 
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Figure 5: Vortex pattern in the fast rotation regime. Pictures of the rubidium condensate 
rotating in the quadrat ic+ quart ic potential, for various stirring frequencies Q/2n. For 
these data oj^jlix = 65 Hz (pictures from [HI])- 

4.4 Transverse monopole mode 

The study of the normal modes of a Bose- Einstein condensate generally provides insightful 
information about the system. In order to gain some understanding of the fast rotation 
regime, we have studied the transverse monopole (or breathing) mode for various rotation 
frequencies of the condensate. 

For a 2D gas at rest in an isotropic harmonic potential of frequency to, this mode 
has a frequency u mp = 2u, which does not depend on the strength of the interactions 
jH21 EH]- The state of the condensate at time t can be derived from the state at time 
by a simple scaling transform. The same result holds for a 3D gas confined in an 
axisymmetric, cigar shaped potential [HI]- In particular the frequency of the transverse 
breathing mode is still u mp = 2uj_ [95]. For a rotating condensate one could have naively 
expected that the frequency of the mode is changed to 2yu^ — f2 2 as a consequence of 
the centrifugal potential. As shown in jnSj this result is not correct and the predicted 
frequency is still u mp = 2u± for all rotation frequencies Q. This is a striking example 
of the influence of the Coriolis force on the system: even though it affects neither the 
BEC transition temperature nor the cloud's equilibrium shape at T = 0, it has a strong 
impact on the condensate's normal modes. Note that the result u mp = 2u± holds for any 
2-dimensional gas with contact interactions confined in a harmonic potential [921 EHl > so 
that the monopole mode cannot be used to monitor any phase transition - like that to a 
strongly correlated state. 

We have studied experimentally this mode in the quadratic+quartic potential de- 
scribed above j9H|- We have checked that the frequency u mp remains at 2u± for rotation 
frequencies notably smaller than u>j_, for which the quartic term plays no significant role. 
When the rotation frequency approaches uj± we measure however a small deviation from 
this value. This deviation increases with the rotation frequency and reaches ~ 10% when 
Q ~ uij_. This result is a consequence of the action of the quartic potential and it can be 
accounted for by a simple analytic model [96 . Due to this deviation the monopole fre- 
quency might represent a sensitive tool to monitor the emergence of new quantum phases 
of the rotating gas. 

A remarkable feature of the transverse monopole mode in the region Q ~ uj± is the 
time evolution of the density profile of the cloud. Instead of being simply a scaling 
transform as in the pure harmonic case, we observe a phenomenon of entering waves 
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Figure 6: Transverse breathing mode of a condensate rotating at — 1.05 ujj_ (Q/2tt = 
68 Hz). The condensate is confined in the quadratic+quartic potential described in the 
text. The structure of the mode corresponds to entering waves, instead of a simple scaling 
transform as in the pure harmonic case. The time interval between two successive pictures 
is 1 ms. 

(Fig. EJ). This structure can be explained by noticing that for Q ~ u)j_ several m = 
modes have a frequency close to u mp [HE], so that beating between them can lead to the 
observed phenomenon. 

5 Conclusions and perspectives 

To conclude, the physics of a rotating Bose gas presents strong analogies with several 
aspects of condensed matter physics: superconductivity in large magnetic fields, Quantum 
Hall phenomena, superfluidity and rotating bucket experiments. It has already led to 
spectacular findings such as the possibility to directly visualize the vortices and to observe 
their vibration modes such as the Kelvin mode (oscillation of a vortex line) and the 
Tkachenko mode (oscillation of the vortex lattice). However, important aspects of the 
problem remain experimentally unexplored. Let us briefly outline three lines of research 
that seem very promising: 

• The possibility to generate quadratic + quartic potential opens the way to the nu- 
cleation of stable giant vortices. In a quadratic potential a vortex with a topological 
charge larger than 1 is unstable, and it can only be observed in a transient way 
[2H1ISZ1- When a quartic potential is present, this instability may disappear and a 
giant vortex can be stabilized at the center of the trap, possibly surrounded with 
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singly-charged vortices [H31 EH EH E1B E3 EH1 EH EU 

• The combination of rotation and optical lattices opens a very interesting class of 
problems. If a ID optical lattice is applied along the axis of rotation, one obtains a 
stack of rotating parallel layers. The structure of vortices in this system remains to 
be investigated. One can qualitatively expect that the vortex cores in neighboring 
layers will remain aligned if the tunneling between layers is large enough, whereas 
they may decorrelate otherwise [HE]- 

• When the rotation speed increases, the number of vortices N v also increases. When 
it becomes of the order of the particle number N, one expects that the ground state 
of the gas will no longer be well described by a mean-field approximation. Instead it 
becomes a strongly correlated state, with a structure very similar to those appearing 
in fractional quantum Hall effect [T31 EH EES EH E] • In practice these type of states 
are expected to be observable only for small particle numbers. 

With such general lines of research still fully open, the next few years should bring us 
a lot of novel and fascinating physics. 
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